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Abstract. The main result of this paper is a bound on the distance between the distribution of an 
eigenfunction of the Laplacian on a compact Riemannian manifold and the Gaussian distribution. 
If X is a random point on a manifold M and / is an eigenfunction of the Laplacian with L'^-norm 
one and eigenvalue —ji, then 



dTvU{X),Z) < -E|||V/(X)||2 -E||V/(X)| 



, This result is applied to construct specific examples of spherical harmonics of arbitrary (odd) 

degree which are close to Gaussian in distribution. A second application is given to random linear 
combinations of eigenfunctions on fiat tori. 
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1. Introduction 



A question which has arisen and been studied often throughout the history of probabihty has 
been to understand the value distributions of certain naturally occurring functions. The simplest 
example of such a question might be to understand the distribution of the number of times out 
, of N trials that a fair coin lands on heads. A more geometric example (studied by Maxwell and 

j. I Borel among others) is understanding the value distribution of the function / : S"~^ M defined 

■ by /(x) = xi. Value distribution in this context refers to endowing the sphere §"~^ with a natural 

probability measure (normalized surface area), and considering the distribution of the random 
lO ' variable which is the push- forward of that natural measure on the sphere to M by the function /; 

. i.e., if X is a random point on the sphere, the value distribution of / is the distribution of the 

^ I random variable f{X) on M. 

A particular class of functions of interest are eigenfunctions of the Laplacian on Riemannian 
manifolds. The main result of this paper is to give a sufficient condition for the approximate 
^ ' normality of these functions, as follows. 

H : 

^ - Theorem 1. Let M be a compact Riemannian manifold (without boundary), and f an eigenfunction 

for the Laplacian on M with eigenvalue —A < 0, normalized so that ^^i^j^-^ Jj^ = 1- Let X be a 
random (i.e., distributed according to normalized volume measure) point of M. Then 

dTvifiX),Z) < ^e|||V/(A)|P-E||V/(A)|P 

where dxv denotes the total variation distance between random variables and Z is a standard 
Gaussian random variable on M. 

That is, the value distribution of / is close to Gaussian if the expected distance between the 
random variable ||V/(A)|p and its mean is small, relative to the size of the eigenvalue. It should 
be pointed out that E||V/(A)|p = A (this is an easy consequence of Stokes' theorem), thus the 
bound on the right-hand side is not always non-trivial. However, it will be seen in examples below 
that there is sometimes sufficient cancellation to obtain a normal limit. 

Eigenfunctions of the Laplacian have arisen as important and natural objects in several branches 
of analysis; in particular, they of central interest in quantum chaos, whose basic problem is the 
understanding of the quantization of a classical Hamiltonian system with chaotic dynamics. See 
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|12| for a readable introduction. A key question in the chaotic case in whether individual eigen- 
functions behave like random waves (see [2]). It is noted in [12] that, were this the case, it would 
imply that the distribution of high eigenfunctions would be approximately Gaussian. In connection 
with this conjecture, Hejhal and various coauthors (see, e.g., [7]) and Aurich and Steiner [Ij have 
provided numerical evidence indicating that on certain hyperbolic manifolds, the distributions of 
eigenfunctions become close to Gaussian as the eigenvalue tends to infinity. Theorem [1] above pro- 
vides a different perspective on the approximate normality of eigenfunctions. In the applications 
given below, asymptotic results are obtained as the dimension of the manifold in question tends to 
infinity. However, there is no a priori reason that Theorem [1] could not also be applied in the case 
of a fixed manifold, as the eigenvalue tends to infinity. 

The results of this paper can be seen as an abstraction of earlier work of the author ^ and 
joint work with M. Meckes [T0|. In those papers, it was shown that the value distributions of 
linear functions on the orthogonal group and on certain convex bodies are close to Gaussian, using 
similar methods to those used here. The proof of Theorem [1] makes clear that eigenfunctions of 
the Laplacian are the natural candidates (at least for the techniques used here) to replace linear 
functions on more general Riemannian manifolds. The main result of [9] on the orthogonal group 
and the results of [lOj in the case of the sphere can be recovered as applications of Theorem [TJ 

The contents of this paper are as follows. Notation and conventions are given in Section 11.11 
below. Section [2] gives the proof of the main theorem, in fact proving a slightly more general 
version which treats linear combinations of eigenfunctions whose eigenvalues are close together. 
Section [3] gives applications of Theorem [1] to the distribution of certain eigenfunctions on the 
sphere and section H] contains applications of the generalization of Theorem [1] from section [2] 

to random eigenfunctions on the torus T*^. The unifying idea in these two cases is that in situations 
in which eigenspaces have high dimension, or the closely related situation in which there are many 
low dimensional eigenspaces whose eigenvalues are very close together, Gaussian distributions occur 
when taking linear combinations of many orthogonal eigenfunctions. 

In the case of the sphere, we have the following result for the second-order eigenfunctions. 

Theorem 2. Let g{x) = j o-ijXiXj^ where A = (aij)fj^i is a symmetric matrix with Tr (A) = 0. 
Let f = Cg, where C is chosen such that \\f\\2 = 1 when f is considered as a function on S"""^. 
Let W = f{X), where X is a random point on S*""^. If d is the vector in M" whose entries are the 
eigenvalues of A, then 

dTv(W,Z) < 

where \\d\\p = ■ 

This theorem gives a fairly complete picture of when second-order eigenfunctions of the spherical 
Laplacian have Gaussian value distributions. The bound on the right-hand side is of the order 
in the case that all of the di are of roughly the same size; e.g., when 



di 



1 i < § 
-1 <i<n 



and n is even. In the opposite extreme case, e.g. for the function /(x) = c{x1 — x\) for a suitable 
normalization constant c, the bound is of order 1. This is what is expected in this case; since it is 
known that xi and X2 are asymptotically distributed as independent Gaussians (as the dimension 
tends to infinity), the random variable f[X) is approximately distributed as Z^ — Zl for independent 
normal random variables Zi and Z2 and X a random point of a high-dimensional sphere. 
In the higher degree case, we have the following theorem. 
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Theorem 3. There is a set o/ orthonormal eigenf unctions of degree £ on ^ and a 

constant c depending on I such that for p : S"^"*^ M defined by 

n 

^aipf{x) 

i=l 

with X] = 1; 

dTvipiX),Z) <cMl 

If the vector of coefficients {aj} is chosen to be random and uniformly distribtuted on the sphere, 
then there is another constant d such that 



mTv{p{x),z) < 



c' 



Si) 



The functions p^ are given exphcity in Section [3l The second statement of the theorem says 

that the typical hnear combination of the p^*^ has an approximately Gaussian value distribution. 
As in the second-degree case, the explicit form of the bound in the first statement of the theorem 
means that for specific choices of a whose components all have roughly the same magnitude; e.g., 
a,- = -1= for each i, the bound is of order -A=. 

In the case of the torus, we have the following result on the distributions of random linear 
combinations of eigenfunctions. 

Theorem 4. Let B be a symmetric, positive definite matrix, and consider the flat torus R^/Z" 
with the metric given by {x,y)B = {Bx,y) . Define a random function f on (¥",5) as follows. Let 
{at,}t,gv be a random point of the sphere of radius \pl in R^, where V is a finite set of vectors of 
M". Suppose that the set V is such that ifv£V, then Bv G Z". Define f by 



fix) = R [ Ea^e^"'^'^'''''M 



Then for fiy = (2-k\\v\\b)'^ and /x > 0, 



EdTv{f{X),Z) < 



8(27r) 



E {v, w)l + (2V2 + V^) ^(^, - /,)2 



Section H] contains a discussion of the bounds above and concrete examples in which they are 
small. The simplest possible case is that of the usual square-lattice torus, in which the matrix 
B = I. In that case, consider the set V to be simply the standard basis vectors of R". Then the 
fact that / as constructed above has a value distribution which is close to normal is essentially just 
the classical central limit theorem. In that case, the second term of the bound above is zero (take 
fi = (27r)^), and the first term is 

In the case of more general B and V, the first term is a weakening of the orthogonality in the 
previous case; it is small as long as V is large but the elements of V have directions which are fairly 
evenly distributed on the sphere. The second term is a weakening of the fact that, in the previous 
example, all of the summands in the linear combination were from the same eigenspace. By taking 

to be a slight perturbation of /, the eigenspaces may become low- or even one-dimensional. 
However, the corresponding eigenvalues may still be clustered very close together, in which case 
the second term of the error above is small. 
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1.1. Notation and Conventions. For random variables X and Y , the total variation distance 
between X and Y is defined by 

dTv{XX) = sup|P(X A)-¥{Y ^ A)l 

A 

where the supremum is taken over measureable sets A. This is equivalent to 

dTv{X,Y) = isup|E/(X)-E/(y)|, 
^ f 

where the supremum is taken over functions / which are continuous and bounded by one. 

Let (M, g) be an n-dimensional Riemannian manifold. Integration with respect to the normalized 

volume measure is denoted dvol, thus /^^ Idvol = 1. For coordinates |^|. on M, define 

\ , g{x) = det(G(x)), g'^ (x) = (G-i(x)),,-. 

X I 

Define the gradient V/ of / : M ^ M and the Laplacian A^/ of / by 

Let <I>*(x,t;) denote the geodesic flow on the tangent bundle TM and let vr : TM ^ M be the 
projection map of a tangent vector onto its base point. 

Acknowledgements. I would like to thank Pierre Albin, Persi Diaconis, Geir Helleloid, Zeev 
Rudnick, Laurent Saloff-Coste, and Steve Zelditch for helpful discussions. In particular, I learned 
about hypergeometric series from Geir Helleloid and about Liouville measure from Steve Zelditch. 



d_ 

dxi 



_d 

dx 



2. The main result 

As mentioned above, the main result of this section is a slight generalization of Theorem [T] of the 
previous section. The idea is that for the method of proof used here, the functions considered need 
only be approximate eigenfunctions of A^; in particular, linear combinations of eigenfunctions all 
of whose eigenvalues cluster about one value are close enough. 

Theorem 5. Let {M,g) be a compact Riemannian manifold without boundary, and let 

be a sequence of mutually orthogonal eigenfunctions of the Laplacian Ag on M with eigenvalues 

—fii < 0. Assume that ffdvol = 1 for each i. Define the function / : M ^ M by 

m 

fix) = "^aifiix), 

i=l 

where the ai are scalars such that = 1. Let X be a random (i.e., distributed according to 

normalized volume measure) point of M . Then for /x > 0, 

dTv(/(X),Z)<- 

where Z is a standard Gaussian random variable on M. 

The idea is to choose /i = — /^j, and consider collections {fi} of eigenfunctions such that the 
eigenvalues jjn are very close together. In particular, for manifolds on which eigenspaces of A have 
high dimension, one can choose all the /ij to be equal, in which case the second error term drops 
out. 



E 



||V/(X)f-E||V/(X)|p 
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The proof is an application of the following abstract normal approximation theorem. A proof of 
the theorem was given in [9] in the case that E' = 0; the proof of the theorem below is a trivial 
modification of the earlier case. 

Theorem 6. Suppose that {W,W^) is a family of exchangeable pairs defined on a common prob- 
ability space with KW = and KW"^ = a"^ . Suppose there are random variables E = E{W) and 
E' = E'iW), deterministic functions h and k with 

lim h(e) = lim k(e) = 

and functions a and (3 with 

^a{a~^W)\ < oo, E|/?(cT-^VF)| < oo, 



such that 
(I) 



(2) 
(3) 



[W, - W\W] = -XW + E' + h{e)a{W), 



^ E [{W, - Wf\W] = 2Xa^ + Ea^ + k{e)P{W), 



£2 



^E\W,-W\^ = 0(1), 



e 



where o(l) refers to the limit as e — > with implied constants depending on the distribution ofW. 
Then 



dTv{W,Z) < J 



^e\e'\ +e\e\ 



where Z ~ m{0,a'^). 



Proof of Theorem [5l To apply Theorem [6] above, start by constructing a family of exchangeable 
pairs of points in M parametrized by e. Let X be a random point of M and let e > be smaller 
than the injectivity radius of the exponential map at X. Since M is compact, there is a range of e 
small enough to work at every point. Now, choose a unit vector V £ TxM at random, independent 
of X, and define 

X, = expxieV). 

To see that this defines an exchangeable pair of random points on M, we show that 

/ g{x,Xe)dfi= / g{xe,x)dfi 
JmxM JmxM 

for all integrable g : M x M ^ M, where is the measure defined by the construction of the pair 
{X,X,). 

Let L denote the normalized Liouville measure on SM, which is locally the product of the 
normalized volume measure on M and normalized Lebesgue measure on the unit spheres of TxM. 
The measure L has the property that it is invariant under the geodesic flow v). See [3], section 
5.1 for a construction of L and proofs of its key properties. By construction of the measure ^, 

/ g{x,Xe)d^= I g{x,TT o ^^[x,v))dL{x,v) 
JmxM Jsm 

= / gi-n- o^'^{y,r]),y)dL{y,r]) 
Jsm 



6 



E. MECKES 



by the substitution {y,r]) = ^'^{x,v), since Liouville measure is invariant under the geodesic flow. 
Now, 

as both correspond to going backwards on the same geodesic. It fohows that 



SM 



g{x, x^)dL{x, v) 



SM 



g{Tr{-^'{y,-7j)),y)dL{y,rj) 
g{-K{-^'{y,r])),y)dL{y,r]) 



SM 



g{TT{'^'{y,r])),y)dLiy,r]) 



SM 



g{ye,y)dfi, 

MxM 

where the second line follows from the substitution —r]^r], under which Lebesgue measure on the 
tangent space is invariant, and the third line follows because Tr{—v) = 7r(u) (both have the same 
base point). Thus the pair {X,X^) is exchangeable as required. 

Now, let f{x) = Y2'^ifii^) ^ statement of the theorem. For notational convenience, let 

W = f{X) and = f{X^). Since {X,X^) is exchangeable, (H^, We) is an exchangeable pair as 
well. 

In order to verify the conditions of Theorem [6l first let 7 : [0, e] ^ M be a constant-speed 
geodesic such that 7(0) = X, 7(e) = X^, and 7'(0) = V. Then applying Taylor's theorem on M to 
the function / o 7 yields 



/(X,)-/(X)=e- 



d(/o7) 



(1) 



dt 



+ 



d2(/o7) 



t=0 



dt2 



+ 0(e3) 



t=0 



e-dxf{V) + - 



d2(/o7) 



dt^ 



i=0 



where the coefficient implicit in the O(e^) depends on / and 7 and dxf denotes the differential of 
/ at X. Recall that dxf{v) = (V/(x),f) for v £ T^M and the gradient V/(x) defined as above. 
Now, for X fixed, V is distributed according to normalized Lebesgue measure on SxM and dxf is 
a linear functional on TxM. It follows that 

E [dxf{V)\X] = E [dxf{-V)\X] = -E [dxf{V)\X] , 
E[dxf{V)\X] =0. 



so 



This implies that 



lim-E[/(X,)-/(X)|X] 



exists and is finite. Indeed, it is well-known that 



limlE[/(X,) - fiX)\X] = ^AJiX) 



for n = dim{M), thus 



lim 4tE \We - W\W] = -— E 



J2(^^^MfiiX)\W 



i=l 



-—W + —E 
2n 2n 



1=1 



W 



It follows that A = ^ and E' = ^E ai(/i - //i)/i(X)| I^] . The higher order terms in e of 
E [Wf — satisfy the integrability requirement of Theorem [6] since / is smooth and M is 
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compact. Finally, 



A 


< 


A* 




< 


1 






m\ 















E 



^ai(/ii - fi)fi{X) 



since the are mutually orthogonal. 
Consider next 

limirE\(W,-Wf\X]. 

By the expansion ([T|), 

E[{W,-Wf\X] = E[{f{X,)- f{X)f\X] 

= e^E[{dxf{V)f\X]+0{e^). 

Choose coordinates |^| ^ neighborhood of X which are orthonormal at X. Then 



thus 



[d.f{v)f = [{Vf,v)f 



V-/ dfV 2 ^ df df 



i=l 



Since V is uniformly distributed on a Euclidean sphere, E[Vilj] = ^Sij. Making use of this fact 
yields 

hm Ie [{dxf{V)f\X] = -\\Vff = ^ + i [||V/f - ^^] , 



thus condition (2) is satisfied with 



E = -E[\\Vf{X)f-^,\W]. 



By Stokes' theorem, 

E||V/(X)f = -E[/(X)A,/(X)] 



E 



f{X)Y,^if^ifi{X) 



jMf{X)+E 



fiX)^a,ifi,-f,)fi{X) 



/i + E 



/(X)5^a,(^,-/i)/,(X) 
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Thus 



n 

_ 1 
n 



l|V/f -/i-E 



f{X)Y,Ci^{^l^-^Ji)f^{X) 



+ E 



f{X)Y,Ci^{^X^-^i)fi{X) 



w 



E[||V/f-E||V/f|VF] +E 



f{x)Y,a^{^^^-^^)fi{x) 



and so 



\^E\ < - 
A 



< 



1^ 



2 

< - 



E 



E 



E 



||V/(X)f -E||V/(X)f 
||V/(X)f -E||V/(X)f 

||V/(X)f -E||V/(X)f 



+ E /(X)^a,(^,-^)/,(X) 

+ v/e72(x) 



E|^^a,(;u,-//)/,(X)^ 



Finally, ([T]) gives immediately that 

E [iVFe - W\^\W] = O(e^). 



□ 



3. The sphere 



Eigenvalues and eigenfunctions of the Laplacian on the sphere S"~^ are well-studied objects. In 
particular, it is discussed in section 9.5 of [H] that the eigenfunctions of the spherical Laplacian 
are exactly the restrictions of homogeneous harmonic polynomials on to the sphere, and that 
such polynomials of degree / have eigenvalue —l{l + n — 2). It follows, then, that the simplest class 
of functions on the sphere to which one could try to apply Theorem [5] are the linear functions. As 
was mentioned in the introduction, understanding the value distribution of linear functions on the 
sphere is an old problem. A careful history is given in the paper of Diaconis and Freedman [5], in 
which the following theorem (and more general versions) are proved. 

Theorem 7 (Diaconis-Freedman). Let X he a random point on the sphere ^/n'S>"^^^ C R". Then 

dTv{Xi,Z) < 

n — 1 

where Z~ 91(0,1). 

It is straightforward to reproduce this result as a consequence of Theorem [5j In Section 13.11 
below. Theorem [5] is applied to the next obvious example, the second order spherical harmonics. 
Section 13.21 treats more specialized examples of spherical harmonics of arbitrary (odd) degree. 

The following lemma, taken from [6], for integrating polynomials over the sphere will be useful 
in applications. 

Lemma 8. Let P{x) = |xi|"i|x2|"^ • • • |x„|"". Then if X is uniformly distributed on S'^~^ , 

r(/3i)---r(/3„)r(§) 



E[P{X)] 



r(/3i + --- + /3„)W2' 



where Pi = ^{ai + 1) for 1 < i < n and 

m = / 

Jo 



s^-^e-'ds 



-"■'dr. 
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3.1. Second order spherical harmonics. 

Theorem 9. Let g[x) = X^i j o-ij^i^j, where A = {aij)^^^^ is a symmetric matrix with Tr (A) = 0. 
Let f = Cg, where C is chosen such that \\f\\2 = 1 when f is considered as a function on S^~^ . 
Let W = f{X), where X is a random point on S^^^ . If d is the vector in whose entries are the 
eigenvalues of A, then 

dTv{W,Z)<VQ^^^'^^^^^ 

where \\d\\p = (Z^.m'^" . 



Proof. To apply Theorem [5l first note that / is indeed an eigenfunction of the Laplacian on 
/S'"~^;since A is traceless, g is liarmonic on and thus is an eigenfunction with eigenvalue A = — 2n. 

Next, observe that g{x) = {x,Ax) , and j4 is a symmetric matrix, so A = U*DU for a diagonal 
matrix D and an orthogonal matrix U. Thus 

g{x) = {x, U*DUx) = {Ux, DUx) . 

Since the uniform distribution on the sphere is invariant under the action of the orthogonal group, 
this observation means that it suffices to prove the theorem in the case that A = diag{ai, . . . , an)- 
Since / is an eigenfunction, the second error term of Theorem [5] vanishes. By the theorem and 
an application of Holder's inequality, 

dTv{W,Z)<l-^^Var{\\Vf\\^), 

and the main task of the proof is to estimate the right-hand side. A consequence of Stokes' theorem 
is that 

/ WVgfda = \X\ I g^da, 

SO one can calculate EHV^p in order to determine C. The spherical gradient \7g is just the 
projection onto the hyperplane orthogonal to the radial direction of the usual gradient. Letting 

stand for summing over distinct indexes, 



EllVc 



E\\VR^g{x)f -Eix-VRug{x)y 



E 



.i=l 



2 2 



4E 



E2 4 I V^' 2 2 



i=l 



n 



i=l 



A "■ 

4 ^-^ . 
— > a: 



n 



i=l 




n{n + 2 



n + 2 



E 

1=1 



ah 



where Lemma [H] has been used to get the third line, and the fact that Tr (A) = is used to get the 
last line. From this computation it follows that the constant C in the statement of the theorem 

should be taken to be 

Now, 



n{n+2) 
2||a||i^ 



, where a = (aj)"^;^. 



E\\Vsn-igf <E\\VRug\\\ 
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since Vgn-ig is a projection of Vigng. 



EllVRnffll 



16 E 



2 2 2 



16 



n{n + 2) 



E4 4 I V^' 2 2 

id 



and so 



4 a 



4n(n + 2) 



4n(n + 2) 



1 + 



2 ||a||| 

kill 



This gives that 



thus 



Var (||V/f) < 8n^ 



+ 8n 



1 + 



2||a| 



n/ V 2 / n V n Wlalh 



since ||a||4 



> n-V4||a| 



2- 



Remark: For n even, consider the quadratic function 

m-- 



n + 2 



j=i 



□ 



This / satisfies the conditions of the theorem and is normahzed such that E/^(X) 
random point on S""^. Theorem [9] gives that 

'6 



1 for X a 



dTvifiX),Z) < a/- 
V n 

so in this case, Theorem gives a rate of convergence to normal. 

The cases in which the bound does not tend to zero as n — > oo are those in which a smaU 
number of coordinates of X control the value of f{X); in those situation one would not expect / 
to be normally distributed. For example, if f{x) = cx\ — cx\, where c is the proper normalization 
constant, then because xi and X2 are asymptotically independent Gaussian random variables (see 
[5]), / is asymptotically distributed as c{Zl — Z|) for independent Gaussians Zi and Z2. 

3.2. Higher degree spherical harmonics. Fix an odd integer i. A canonical example of an 
eigenfunction of the spherical Laplacian of degree (. is the polynomial pi{x) given by 



Pi{x) = AC^ 2 (x„), 
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where A"^ = ^"^^^^'Ig";^! 2)^ ^ normalization constant chosen such that Kp'j(X) = 1 for X a 
random point of S^~^, and 

is a Gegenbauer polynomial. The eigenvalue of this eigenfunction of the Laplacian on is 
i{n + e-2). 

The reason that this particular eigenfunction is natural is that is arises in the representation 
theory of the orthogonal group. Let T"^ be the irreducible representation of SO{n) in the space 
<^n/ 1 ^2i>>p^n,i-2 ^ where ^H*^'^ is the space of homogeneous polynomials of degree I van variables. The 
polynomial pi{x) above is the zonal polynomial of degree I of the representation T"^ relative to the 
subgroup SO{n — 1), where as usual, SO{n — 1) is identified with the subgroup of SO{n) fixing the 
north pole of S'^~^. See [14] for background on spherical harmonics and Gegenbauer polynomials, 
including all facts stated above. 

/ \ n — 2 /. \ 71 — 2 

Observe that since {x) = AC^ ^ is an eigenfunction, so are p^ (x) = AC^ ^ (xk) for all 
1 < k < n. Let ai, . . . , On C R such that X] = 1; ^^'^ consider the degree £ eigenfunction 



k=l 

n — 2 

Note that Ep^(X) = 1 for X a random point of S"""^, because if £ is odd, then the ^ (x^) are 
orthogonal for different k. In this example, a bound on the total variation distance between the 
distribution of the random variable p{X) for X random on the sphere and a standard Gaussian 
random variable Z is obtained in terms of the Oj. Specifically, 

Theorem 10. There is a constant c, depending on £, such that for p : S"~^ — > M defined by 

n 

^aip''f{x) 

i=l 

with — ^! 

dTvipiX),Z) <cMl 

If the vector of coefficients {aj} is chosen to be random and uniformly distribtuted on the sphere, 
then there is another constant d such that 

mTv(j>{x),z) < ^. 



Remark: The second statement of the theorem says that, typically, spherical harmonics ob- 
tained as linear combinations of the p^^\x) are approximately Gaussian in distribution. As in 
Theorem [9l it may also be helpful to consider a specific example. If |aj| = for each i, then 



dTv{p{X),Z) < 



for a constant c depending only on £. As discussed after the statment of Theorem[9l the bound ||a||4 
is not small in the case that only a fixed small number k of Oj are non-zero. In such a case, p{X) 
is not distributed as a Gaussian random variable but as a degree £ polynomial in k independent 
Gaussian random variables. 
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Proof. Recall that Theorem [5] gives 



1 



(3) 



since 



dTvip{X),Z) < -VVar(||V5n-ip(X)||2) 



sn-ip{X)f =i{n + i-2) 
by Stokes' theorem. Thus to apply the theorem, we must compute (or bound) 

Consider the function p{x) = X]fc=i ^kpf^\x) on a neighborhood of S"'^^. Then for x G S"'~^, 

\\Vsn-ip{x)\\' = ||Vp(x)f - {x,Vp{x)f , 

where Vp{x) denotes the usual R"-gradient of p. Thus 

E||Vsn-ip(X)f <E||Vp(X)|h 

bounding the right hand side turns out to be enough. 

It is proved in [14J that the Gegenbauer polynomials have the property that 

d 



(4) 
thus 



dt 



C^(i) = 2pC^+_\(i), 



Vp{x) = [Aatin - 2)C^^_^^{xk) 
It is therefore necessary to estimate 



k=l 



E 



(5) 



fc=i 



n 

Y,4iclA^k))'+Y^ala]{cl,{x,))HcL{xj)f 



A^(n-2rE 



_k=l k,j=l 

where stands for summing over distinct indices. 
The proof hinges on the following key facts. 

(1) For k + j, 

A\n - 2)^E [{Cl,{xk))\Cl,{x,)f\ =f{n + 0(1))'. 

(2) There is a constant c such that for all k, 

Using the normalization X] ^fc ~ these facts imply that 

n 

E||Vsn-ip(X)f < £V + cV J^^at + c" 



n 



k=l 



for constants c/, c". Note that since ||a||2 = 1, ^cif. > ^ and so the last term can be absorbed into 
the middle term. The result is now immediate from ([3|). 
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E 



It remains to verify the two key facts above. By ([2]), 

.2 

4 



r(i) 



E 

k,m,p,q=0 



Lemma [8] gives 



.i)k+ni+p+qY (£_i + |_A;)r(£-i + §-m)r(^-i + §-p)r(^-i + §-g) 

22fc+2m+2p+2g;.!^!p!g!(£ - 1 - 2/c) ! - 1 - 2m) ! (£ - 1 - 2p) ! - 1 - 2^) ! 

E 



jl^2£-2fc-2m-2ji^2£-2p-2q-2 



E 



2^-2fc-2m-2j^2^-2p-2<j-2 



r(«-fc-m-i)r(<-p-9-i)r(2) 



ir 

vr 



vrr (f + 2£ - - m - p - g - 2) 



1 



k-m--]T{e-p-q--]{- + 0{^ 



n 



1 



2l-k-m-p-q-2 



where here and throughout, O (^) , 0(1), etc. refer to the hmit as n ^ oo, with ^ still fixed and 
the constants implicit in the O depending on H.. It follows that 

2 



E 



C/_,(Xi)C/_i(X2) 



£-1 
i2£-2 2 



(6) 



^ (-l)^+'"r(^-l + §-A;)r(^-l + f -m)r(£-A; 



fc,m=0 



22fc+2m^!^!(^ - 1 - 2k)\{l - 1 - 2m)! 



n 



'k—m—l^ 



Now, 
(7) 



n 



r^-i + 7T-^=rU 0+0(1) 



n\ /n 



e-k-i 



so the right-hand side of ([6|) reduces to 



2 



_^ ^_l)fe+r«r(^- A: 

(2n + 0(l)) 22fc+2™A:!m!(£ - 1 - 2k)\{e - 1 - 2m)! 



m-i)7r-V2 



fc,m=0 



Recall that 



thus 



yl4(n-2)% C/„i(Xi)C/„i(X2 



2 



(^!)^(n + 0(l))M2^-i ^ 



(-i)^+™r {i-k 



m-i)7r-V2 



fc,m=0 



22fc+2m/.!^!(£ _ 1 _ 2A;)!(£ - 1 - 2m)! 
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Claim: For £ G N odd, 



(-i)'=+™r {i-k 



m 



-1/2 



1 



22fc+2mA;!m!(^ - 1 - 2A:)!(£ - 1 - 2m)! - 1)! ' 



k,m=0 

A proof of this claim is given in the appendix. Making use of the claim, it is immediate that 

r — H "I 2 ri 

A^n - 2)4E [c/_,(Xi)C/_,(X2)J =e{n + 0(1))'. 
It remains to verify fact 2: there is a constant c such that for all k, 



< cn 



In fact, this is much easier than the previous part since the exact asymptotic constant is not needed. 
Proceeding as before, 



E 



C^/-i(^i) 



f-i 



E 

k,m.,p,q=0 



22fc+2'»+2p+2g/^!^!p!g!(^ - 1 - 2fc)!(^ - 1 - 2m)!(^ - 1 - 2p)!(^ - 1 - 2g)! 

E 



^1 



From Lemma [HI 



E 



4^-2fc-2m-2p-25-4 
1 



r (2£ - A; - m - p - Q - I) r (f ) 
vrF (§ + 2£ - fe - m - p - g - 2) 

-r ( 2^ - A: - m - p - g - - ) ( - + O ( 4 



2 / V n 



and so 



E 



C/_i(^i, 



2£-fc-m-p-g-2 



-,4£-4 



E 



A:,r?i,p,(j'=0 

making use again of d?]). Since 



(_;L)fc+m+p+gp (2£-A;-m-p-g-|)(f+0 (1)) ' 
22fc+2m+2p+2q^!^!p!g!(^ - 1 - 2A;) ! - 1 - 2m) ! - 1 - 2p) ! - 1 - 2g)!7r I ' 



= (^!)2(n + 0(l)) 



-2^ 



it follows that 



n 



where C (which depends on £) is independent of n. 



□ 
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4. The torus 

This section deals with the value distributions of eigenfunctions on fiat tori. The class of functions 
considered here are random functions; that is, they are linear combinations of eigenfunctions with 
random coefficients. 

Let {M,g) be the torus T" = M^/Z*^, with the metric given by the symmetric positive-definite 
bilinear form B: 

{x,y)B = {Bx,y) . 
With this metric, the Laplacian on T" is given by 

A«/(.)=j:(is-'),,^(.). 

Eigenfunctions of Ab are given by the real and imaginary parts of functions of the form 

for vectors v G M" such that Bv has integer components, with corresponding eigenvalue — /Li^ = 
Consider a linear combination of eigenfunctions with random coefficients: 

fix) :=3fi(^a,e2-<^"'-)|, 
\vev 



where V is a finite collection of vectors v such that Bv has integer components for each v eV and 
{0'v}veV is a random vector on the sphere of radius ^/2 in M^. Assume further that v + w ^ for 
v.w € V. Then the following theorem holds. 

Theorem 11. For a random function f on (T", B) defined as above and for fJ, > 0, 

1 



mTvif{X),Z) < 



8(27r)^ 



The idea is that the eigenvalues —fiy = —{2Tr\\v\\B)'^ corresponding to the summands be clustered 
about the value fi and that the size of V be large. In some cases, e.g. for B = I, it may be possible 
to eliminate the second error term by choosing a large set V such that all elements have the same B- 
norm, that is, it may be that the eigenspaces of the Laplacian A^ have large dimension. However, 
this need not be true in general. The presence of the second error term is to take into account that 
while it may not be the case that eigenspaces are large, it may be that there are many eigenspaces 
of low dimension whose eigenvalues are all close together. 

Note that requiring the first term of the error to be small is a weakening of an orthogonality 
requirement on the elements of V. By the Cauchy-Schwarz inequality. 



iv,w)'B< < (max||v||^ ) |V| 



v,wGV v,wGV 

If all of the vectors v ^ V have /i^, ~ ^u, then this last expression is approximately (2^)4 1"^^; 

the first error term is bounded by 2\/2. This is of course worthless, since total variation distance 
is always bounded by one. On the other hand, if the elements of V are mutually orthogonal, then 

v,wev vev ^ ' 
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and then the first error term is bounded by -^^=, which is small as long as V is large. The idea is 

to aim for a set V somewhere in between these two extremes: not to require the elements of V to 
be orthogonal, but to choose them to be fairly evenly distributed on the i3-norm sphere on which 
they lie. 
Examples. 

(1) Let B = I. Choose V to be the set of vectors in R" with two entries equal to one and all 
others zero. Then l^l = (g) ^-nd \\v\^ = 2 for each v ^ V. Given a fixed f G V, there are 
exactly 2n — 3 elements u; G V (including v itself) such that (f , w) ^ 0. Thus 

1 v-^ , v9 1 f'n\,.., , 16n — 12 



1^1(1^1 + 2) „i:lv ' ' 0(0+2) \2, 

thus in this case there is a constant c such that 

mTv{f{x),z) < 

\/n 



(2) Let B = diag{\ + (5i, . . . , 1 + with 5i small real numbers. Let V be the set of vectors 
of the previous example. For u G V, let v' = ^jq^vi, . . . , j^^Vr^; v' is constructed such 
that Bv' has integer components. Let V = {v' : w G V}. 



Now, let e = max,- 



1+5,; 



Then 



1^^ lis - 2 



n 

2 



< e\\vf = 2e. 



=1 ' ' i=i 

By construction, (v,w)b 7^ if and only if {v,w) / 0, and in those cases, 

{v,w)l < \\v\\l\\w\\l < (2 + 2e)^ 

Taking /u = 2(27r)^ in the theorem, it follows that for / defined as above using this matrix 
B and index set V , there is a constant C so that 

mTv{f{X),Z) < Cmax ( ) . 



n 



Proof of Theorem\ll[ To apply Theorem [U an expression for HV^/H^ is needed. 



Vb/(x) = \ K I j;(27ri)a,(Bi;),(i?-i),fce2-<^-'- 

, j=l v€V 




using the fact that B is symmetric. 
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It follows that 



(8) 



I|Vb/(x)||| = B,k (Vb/(x))^. (Vb/(x)), 

j,k \v&v J Vtuev / 



2 



2Tri{Bv-Bw,x) _ ^2iTi{Bv+Bw,x) 



Let X be a randomly distributed point on the torus. The first half of the error term in Theorem [5] 
is 



E 



\\VBf{X)\\l-Ex\\VBf{X)\\l 



<jE\\VBf{X)\\%-E{Ex\\VBf{X)\\iy 



Start by computing E||Vb/(X)||^. From above, 



\\VBfix)\\% = 27r^^ 



ayawa^'a^ui{v,w)Biv',w')B 

.w,v',w' 

^^2TTi{Bv-Bw-Bv' +Bw' ,x) _ ^2m{Bv-Bw-Bv' -Bw' ,x) _|_ ^2m{Bv~Bw+Bv' ~Bw' ,x) 
_ ^2ni{Bv-Bw+Bv'+Bw',x) _ ^2ni{Bv+Bw- Bv' +Bw' ,x) _|_ ^2ni{Bv+Bw- Bv' - Bw' ,x) 



^2ni{Bv+Bw+Bv' -Bw' ,x) _|_ ^2ni{Bv+Bw+Bv' +Bw' ,x) 



Averaging over the coefficients {a^} gives 



^a\\VBfix)\\% 



|V|(|V|+2) 



3f? 



_|_g47rj(BD-_Bi'',x) _|_ ^4TTi{Bv+Bv' ,x) 



_|_ 2 ^2 — 2e^'^^^^'"'^'^ — Q-^'^^i^'^'^) — Q'ini(Bw,x) 

_^^4TTi(Bv— Bw,x) _j_ ^4:'!Ti{Bv+Bw,x)^ 
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Now averaging over a random point X of the torus gives 



mBf{x)rB 



|V|(|V| + 2) 
|V|(|V| + 2) 



v,w£V 



Next, 



E 



Ex||Vb/(X)||| 



E 



(27r 



\4 



|V|(|V|+2) 
|V|(|V|+2) 



I ii2 II ii2 



v,w£V 



-III. 



thus 



E\\VBf{X)\\%-E{Ex\\VBf{X)\\lY 



< 



\ 



|V|(|V|+2) 



2 E -^E 



\B 



< 



2{2ttY 



|V|(|V|+2) 



E (■^'^^ 



2 

B- 



Finally, the second error term of Theorem [5] is a multiple of 



vev 



where iJ,y = (27r||'t;||B)^ and ^ = |^ J2vevf^u- Holder's inequality. 



y vev V '^ev 

This completes the proof. 

5. Appendix 
This section contains the proof that for £ odd. 



(9) 



2-^ E 



k,m=0 



_ ^ _ ^ _ 1) vr"i/2(^ - 1)! _ 
22fe+2m/.!^!(^ _l_2k)l{e-l- 2m)! ~ ^' 



□ 
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For n € N, let nil = n{n — 2)(n — 4) • • • (1, 2), where the last factor is 1 or 2 depending on the parity 
of n. Using the fact that T{x) = (x — l)r(x — 1) together with T = ^/^T, the left-hand side of 
& is 



-1 

— (-l)'=+'"(2£-2fc-2m-3)!!(^- 1)! 
2fc+m^!^!(£ -2k- 1)\(£ -2m- 1)1 

k,m=0 ^ ^ ^ ' 

2 / 1 \k+m 



E 



(10) 



(_l)fe+m(2^ -2k -2m- 2)\{l - 1)! 



k m-o 2''+"fe!m!(£ -2k- 1)!(^ - 2m - 1)!(2£ - 2A; - 2m - 2)!! ' 



Next, let p = ^-!^. Then ([T0|) becomes 

P (_l)fc+m(4p_2A:-2m)!(2p)! 
^^^2^+'"A;!m!(2p - 2k)\{2p - 2m)!(4p - 2A; - 2m)!! 

^2-2p V (-l)fc+"^(4p-2fc-2m)!(2p)! 

fc^o ^'"^K2p - 2A;)!(2p - 2m)!(2p - k - m)\ 

.-2pY- (-l)'(2p)! /y^ (-ir(4p-2fc-2m)! \ 
^ k\{2y - 2k)\ I m!(2p - 2m)!(2p -k-m)\y 

and the goal is to prove that this last expression is equal to 1. We do this by showing 

-l)'^(4p-2/c-2m)! _{2^i' k = 



(11) E 



^^j^m!(2p- 2m)!(2p- /c-m)! lo k = l,...,p. 

This can be done by a computer algebra system such as Maple; the following direct proof is due to 
Geir Helleloid [8j. 

Let {a)n be the shifted factorial 

{a)n = a{a + 1) • • • (a + n — 1). 

The following two straightforward identities involving shifted factorials will be useful: 

(-l)'=n! 



{n-k)\ 



{-n)k 



and 



The hypergeometric series 



is denoted 



(ai)n • • • {ap)n 
(^l)n--- (^g)n n! 



E 



pFq = (ai, ...,ap;bi,... ,bq;x) 
with the X omitted if x = 1. The Chu-Vandermonde identity is 

TP f h \ (C + 

2-Fi(-n,-6;c) = —— . 

\C)n 
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The proof of (jlip is as follows: 

{-l)^{Ap-2k-2my. ^ (-l)"^(-l)2™(4p-2A;)!/(-4p + 2A;)2m 



mU')n — ')m\U')'n — h — m\\ ^-^ 



m\{2p - 2m)\(2p - k - m)\ m!((-l)2™(2p)!/(-2p)2™)((-l)'"(2p - k)\/{-2p + k)m) 



{Ap-2k)\ j^{-2p)2^{-2p + kl 



{2p)\{2p - k)\ m\{-Ap + 2k)2m 

{Ap-2k)\ ^22-(-pV(^) {-2p + k\ 



(2p)!(2p - k)\ ^!22m(_2p + kU (- 
{^p-2k)\ ^(-^)-("^)^ 



-4p+2fc+l 
2 



(2|,)!(2p-fe)!^^^, (^ -4p+2fc+i 

(4p-2fc)! ^ (-j^)m(^ 
(2p)!(2p-/c)!^^^, r-4p+2fc+i 



(4p-2/!;)! ^ / -2p+l -4p + 2A; + l 
2F1 -p, 



{2p)l{2p-k)\ V 2 
(4p-2A:)! 



(2p)!(2p - A;)! / -4p+2fc+i 
(4p-2/c)! (/c-p)---(A:- 1)(-2)P 



(2p)!(2p - A;)! (4p - 2A; - l)(4p - 2/c - 3) • • • {2p -2k + l)' 

Now, for A; = 1, . . . ,p, this last expression clearly vanishes. For k = 0, the last line reduces to 

(4p)!(p)!2f _ (4p)!(p)!2P(2p-l)!! 

{2p)\ (2p)! (4p - l)(4p - 3) • • • (2p + 1) ~ (2p)! (2p)! (4p - 1)!! 

23p(p)!(2p- 1)!! 



(2^,)! 

22P. 



This completes the proof. 
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